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ABSTRACT 

Results are presented of a search for cosmic sources of high energy neutrinos 
with the ANTARES neutrino telescope. The data were collected during 2007 
and 2008 using detector configurations containing between 5 and 12 detection 
lines. The integrated live time of the analyzed data is 304 days. Muon tracks 
are reconstructed using a likelihood-based algorithm. Studies of the detector 
timing indicate a median angular resolution of 0.5 ± 0.1 degrees. The neutrino 
flux sensitivity is 7.5 x 10~^ (£',,/GeV)~^ GeV~^ s"-^ cm~^ for the part of the sky 
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that is always visible [S < —48 degrees), which is better than limits obtained by 
previous experiments. No cosmic neutrino sources have been observed. 

Subject headings: astroparticle physics; cosmic rays; neutrinos 



We dedicate this paper to the memory of our colleague and friend Luciano Moscoso, who 
passed away during the preparation of this paper. 



Introduction 



High energy cosmic rays permeate the universe, but their origins remain uncertain. 
Many types of acceleration sites have been sugges t ed, su ch as supernova remnants, micro- 
quasars and active galactic nuclei (see e.g. iBeckerl ( 120081 ) for a review). If such objects are 
sources of cosmic rays, high energy neu trinos may also be produced. These may be detected 
as a diffuse flux (lAguilar et al.ll2011al ). but the observation of point-like sources of cosmic 
neutrinos would offer a unique opportunity to identify and study the sites and mechanisms 
of cosmic ray acceleration. 

In this paper, results are presented of a search for cosmic point-like sources of high 
energy muon neutrinos performed with the first data taken by the ANTARES telescope. 
Located in the Mediterranean Sea, the ANTARES experiment is predominantly sensitive 
to neutrinos from the southern hemisphere in the TeV to PeV energy range. In particular, 
this allows the study of Galactic sources and complements the IceCube observatory, which 
primarily explores the northern hemisphere. 



2. Data Collection 



The detection principle relies on measuring Cherenkov light emitted by high energy 
charged particles that result from neutrino interactions inside or near the instrumented vol- 
ume. In particular, charged current interactions of muon-neutrinos produce high energy 
muons. The muon direction closely follows that of the neutrino and can often be recon- 
structed with sub-degree accuracy. The large background from downgoing muons due to 
cosmic ray interactions in the atmosphere is reduced by selecting only upgoing muons as 
neutrino candidates. 
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The Cherenkov light is detected by an array of photomuhip lier tubes, each hou sed in 
a pressure resistant glass sphere called an Optical Module (OM) ( lAmram et al.ll2002[ ). The 
OMs are placed in storeys of triplets along vertical detector lines, which are anchored to the 
sea bed at a depth of 2475 meters with a spacing of approximately 60 meters. Each of the 
12 detector lines contains up to 25 storeys, separated by 14.5 meters. 



The arrival time and charge of the photomultiplier signals are digitized (lAguilar et al. 



20101 ) into 'hits' and transmitted to shore. An online filter identifies events containing muons 
within the data stream that otherwise consists of optical backgrounds due to natural radioac- 
tive decays and bioluminescence. The selected events are stored for offline reconstruction. A 
detailed d escription of the dete c tor, the data acquisi tion and the online filtering algorithms 



IS given m 



Ageron et all (120111 ): lAguilar et all ( 120071 ) 



2.1. Data Taking and Processing 

The data presented here were collected between January 31st, 2007 and December 30th 
2008. During this time, the construction of the ANTARES detector was still in progress. 
The detector consisted of 5 lines for most of 2007 and of 9, 10 and 12 detector lines during 
2008. The total live time of the data used for the analysis is 304 days.( 144, 38, 48 and 74 
days with 5, 9, 10 and 12 detector lines respectively). 

The online event selection identifies triplets of OMs that detect multiple photons, either 
as a high-charge hit, or as hits separated by less than 20 ns on adjacent OMs. At least 
4 of such triplets are required throughout the detector, with the relative photon arrival 
times being compatible with the hypothesis of the light being emitted along the track of a 
relativistic particle. 



The arrival times of the hits are calibrated as described in lAguilar et al.l ( 1201 Ibl ). The 
inter-line timing has been measured in an iterative procedure by comparing the expected hit 
times to those measured in a large sample of reconstructed downgoing muons. In addition, 
a similar, independent measurement has been done using light from the optical beacon 
system, with compatible results. The positions and orientations of the OMs vary because of 
the sea currents. An acoustic positioning system, combined with compasses and tiltmeters 
located along the detector lines, measures the positions and orientations of the OMs, with 
an accuracy of ~10 cm. 
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2.2. Event Reconstruction and Selection 



From the timing and position information of the hits, muo n tracks are reconstructed 
using a multi-stage fitting procedure, based on iHeijboerl ( 120041 ) . The initial fitting stages 
provide the hit selection and starting point for the final fit. The final stage consists of a 
maximum likelihood fit of the observed hit times and includes the contribution of optical 
background hits. In order to maximize the robustness against inaccuracies in the simulation 
of the optical background, the charge information of the hits is not used in the present 
analysis. 

Neutrino candidates are selected from upgoing events using criteria that have been 
determined in a 'blind' manner; i.e. before performing the search analysis on the data. 
The criteria are chosen to optimize the sensitivity, which is defined as the expected median 
value of the upper limit on the neutrino flux. The angular uncertainty obtained from the 
muon track fit is required to be smaller than 1 degree. Tracks are also required to have a 
good reconstruction quality. The latter is quantified by a param eter A, which is based on the 
value of the likelihood function obtained for the fitted muon (see IHeijboerl ( 120041 ) for details). 
The cumulative distribution of A for muons which are reconstructed as upgoing is shown 
in Figure [1] along with the simulated con tributions froni atm ospheric muons and neutrinos. 
The simulation uses the 'Bartol' model ( lAgrawal et al.lll996l ) fo r the atmospheric neutrino 
flux. The atmospheric muons are simu lated by the QGSJET (IKalmykov fc Ostapchenko 
19931) and CORSIKA (IHeck et al.l Il998l ) packages with the primary cosmic ray flux from 
Nikolskv et all (Il984l ). 



The flnal sample of neutrino candidates consists of 2190 upgoing events with A > —5.4 
and an angular uncertainty < 1 degree, which were selected out of a total of ~ 100 million 
(down going) reconstructed muons. The simulation indicates that the selected sample con- 
tains about 60% atmospheric neutrinos; the rest are mis-reconstructed atmospheric muons. 
The agreement between data and simulation in the selected s ample is well within the ~30% 
systematic uncertainty on the neutrino flux ( iBarr et al.ll20061 ). 



3. Detector Performance 

The angular resolution and effective area of the detector for neutrinos passing the selec- 
tion cuts have been determined using simulations and are shown in Figure [2l The simulations 
include the stochastic energy loss of muons via electromagnetic shower production, the prop- 
agation of Cherenkov photons through sea water and a detailed simulation of the detector 
electronics. 
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3.1. Angular Resolution 

The cumulative distribution of the angle between the reconstructed muon direction and 
the neutrino direction is shown in Figure [2] (left panel) for neutrino events with a neutrino 
spectrum proportional to E'"^, where E^, is the neutrino energy. The median of this angular 
error is 0.5 ± 0.1 degrees. For the subset of data in which the full 12-line detector was 
operational, the resolution is estimated to be 0.4 ± 0.1 degrees. 

The systematic uncertainty on this quantity has been estimated by varying the time 
resolution of the OMs in the simulation. The allowed range of is determined by 
requiring that the A distribution in the resulting simulation be compatible with the observed 
atmospheric neutrino events. The best agreement between data and simulation is obtained 
for At = 2.5 ns. Hence, this value is used for all simulations in this analysis, in particular for 
extracting the central value of the allowed range of angular resolutions. A time resolution 
of 3.4 ns is found to be incon ipatible with the neutrino flux model at the 2a level, where 



the uncertainty is taken from iBarr et al.l (|2006[ ). This places an upper bound on the time 



resolution, which translates into a 1 a systematic uncertainty on angular resolution of 0.1 
degrees. This uncertainty incorporates, to first order, all effects which have a net result 
of degrading the time resolution, such as possible mis-alignments and inaccuracies in the 
simulation of light propagation in the water or the transit time distribution of the PMT. 
A similar analysis with analogous results has been performed using downgoing muon data 
instead of upgoing neutrino candidates. 

The absolute orientation of the detector is known with an accuracy of about 0.1 degrees 



( iHalladj iani |20 lOf ) : this uncertainty is taken into account as an independent effect. 



3.2. Acceptance 

The effective area for muon neutrinos A*^ is defined as the ratio between the selected 
neutrino event rate and the cosmic neutrino flux. It is determined from simulations and is 
shown in Figure [2] (right panel) as a function of the neutrino energy for three declination 
intervals. Throughout this paper, the cosmic neutrino flux is assumed to consist of an 
equal amount of z/^ and z/^. 

In the search, limits are set on the constant in the flux parametrization dN/dE^, = 
X [E^/GeV]^'^ GeV^^cm~^s~^. The acceptance A for such a flux is defined as the constant 
of proportionality between the number of selected signal events and the flux intensity 0. 
It can be computed by convoluting A'^{E) with dN /dEy{E). For declinations 5 < —48 
degrees, A and A'^J^ are approximately constant. For —48 < 5 < 48 degrees, the functions 
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°-7 -6.5 -6 -5.5 -5 -4.5 -4 -3.5 -3 

quality variable A 



Fig. 1. — Cumulative distribution of the reconstruction quality variable, A, for upgoing 
events in data (black points) which have an error estimate < 1 degree. The Monte Carlo 
simulated contributions from mis-reconstructed atmospheric muons (dashed blue line) and 
atmospheric neutrinos (solid red line) are also shown. The bottom panel shows the ratio 
between data and simulation. The vertical line indicates the analysis cut A > —5.4. 
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decrease because of the requirement that the tracks are upgoing. For a source dechnation 
of -90 (0) degrees, A = 3.2 (1.8) xlO^ GeV cm^ s. This means a total of 3.2(1.8) neutrinos 
would be detected and selected from a point source with a flux of 10~^ GeV~^cm~^s~^. For 
this flux model, the energy of 80% of the selected signal is in the range 3 < -E^ < 700 TeV. 

To constrain the systematic uncertainty on the acceptance, the atmospheric neutrino 
data have been compared to a simulation in which the efficiency of each OM is reduced. The 
maximal reduction for which data and simulation are still compatible corresponds to a 15% 
reduction in the acceptance for signal neutrinos with an E'"^ spectral shape. Hence, a 15% 
systematic uncertainty on the acceptance is used in the limit calculations. 




Fig. 2. — Performance of the detector for the selected signal events as determined from sim- 
ulation. The cumulative distribution of the angle between the reconstructed muon direction 
and the true neutrino direction is shown for selected signal events with an E~'^ spectrum (left 
panel). The neutrino effective area A*^, averaged over three declination ranges, is shown as 
a function of the neutrino energy (right panel). 



4. Search Method 

Two distinct approaches to look for point-like neutrino sources have been used. In 
the full sky search, a search is made for an excess of events over the atmospheric neutrino 
background anywhere in the field of view. In addition, a candidate list search is made, where 
the presence of a signal is tested at the locations of 24 known pre-defined high energy gamma 
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ray sources which could be expected to emit neutrinos. They include supernova remnants, 
microquasars and BL Lac objects. While more restrictive in scope, the candidate list search 
requires less signal to reach a significant excess, compared to the all-sky search. The addition 
of the candidate search may therefore be thought of as a way to give extra a-priori weight 
to the preselected sources. 



4.1. Event Likelihood 

The search method is based on the likelihood of observing the events, which is defined 

as: 

\ogCs+h = ^log[/isig X J^{l3i{5s,as)) + B{6i)] - /isig - A^bg (1) 

i 

where the sum is over the neutrino candidate events, and J-" is a parametrization of the 
point spread function. This is defined as the probability density to find the reconstructed 
muon i an angle (3 away from the declination Sg, and right ascension «<, of the source; it is 
closely related to the angular resolution (see Figure |2]). B{6) is a smooth parametrization of 
the background rate derived from the observed declination distribution of the 2190 selected 
events. The mean number of selected signal events produced by the source is /isig. The 
term A^bg represents the total number of expected background events, which is constant and 
therefore does not infiuence the maximum likelihood fits or the likelihood ratio. 

In the candidate list search, the likelihood is maximized for each candidate by numer- 
ically fitting the source intensity /isig to the events located within 20 degrees of the source, 
with the source coordinates fixed to the known position. In the full sky search, potentially 
significant clusters are first identified using a loose cone selection, which requires at least 
4 events in a cone of 3 degrees diameter. For each cluster, the likelihood is maximized by 
fitting the source coordinates and the intensity, yielding maximum-likelihood estimates for 
these quantities. 

The next step is to compute the test statistic, which is defined as the logarithm of the 
likelihood ratio: 

Q = log/:r+b -log£b, (2) 

where is the maximum value of the likelihood found in the fit and £b is the likelihood 
computed for the background-only hypothesis (/isig = 0). A large (small) value of Q indicates 
that the data are compatible with the signal (background). 
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4.2. Statistical Interpretation 

The statistical interpretation of the search result relies on simulated pseudo-experiments 
in which the background events are randomly generated by sampling the declination from 
the parametrization B and the right ascension from a uniform distribution. Events from a 
neutrino point source are simulated by adding events around the desired coordinates accord- 
ing to the point spread function. The systematic uncertainties on the angular resolution and 
orientation of the detector are incorporated by varying the simulated characteristics of the 
signal events within the assigned uncertainties. 

Distributions of Q are obtained from the simulated experiments by applying the search 
method to each of them. Significances in the form of p-values are derived from the Q 
distribution for background-only pseudo-experiments. The distributions including simu- 
l ated signal events are use d to set flux limits following the Feldman- Cousins prescription 



(IFeldman fc Cousins! Il998l ). The uncertainty in the acceptance A is accounted for in the 



limits by numerically integrating over the allowed range defined by the ±15% uncertainty. 



5. Results 

Figure [3] shows a sky map of the selected events in equatorial coordinates along with 
the locations of sources from the candidate list (circles) and of the most significant cluster 
found in the full-sky search (square). 

In the full-sky search, no significant clusters of neutrino candidates have been found. The 
most signal-like cluster of events has a maximum likelihood estimate of the source position 
{as, Ss) = (43.21°, —0.50°). The fit assigns 3.4 events as signal events and finds a value of the 
test statistic Q = 6.8. Such a value, or larger (more signal like), occurs with a probability 
p = 88% anywhere in the field of view amongst background only pseudo-experiments. The 
result is thus compatible with the background only hypothesis. 

The results of the candidate list search are shown in Table [H None of the candidate 
sources have a significant excess of events at their location. The most signal-like source 
candidate is HESS J1023-575, where 3(5) events are within 1(3) degrees of its position. 
For this cluster of events, Q = 2.5. In the absence of signal such a value or larger, has a 
probability of p = 17% to occur amongst the 24 source candidates. 
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5=+90° 




5=-90° 



Fig. 3. — Map in equatorial coordinates of the 2190 selected neutrino candidates. The 
position of the most significant cluster (see text) is indicated by the square. The circles 
denote the positions of the 24 sources from the candidate list. 
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5.1. Neutrino Flux Upper Limits 

As no significant point sources are observed, 90% confidence level limits are obtained 
for the intensity, 09o%cl^ neutrino flux from each of the source candidates. They 

are listed in Table 1 and are shown in Figure H] as a function of the source declination. 
Figure S] also shows the sensitivity of this analysis, defined as the median value of the 
expected limit. It is in agreement with the median value of the actually observed limits. 
For the area of the sky that is always visible, the {6 < —48 degrees) sensitivity is about 
7.5 X {E^/GeVy^ GeV-^ s'^ cir-^. 
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Fig. 4. — Limits set on the normalization of an E^j'^ spectrum of high energy neutri- 
nos from selected candidates (see Table [T]). The points show the 90% C.L. limit at the 
declination of the candidate source. In addition to the present result, several previously 
published limits on so u rces in both t h e sou t hern and nort h ern h emis phere are shown from 
Ambrosio etaP J200l[ ): iThrane et all koO^ : lAbbasi et all toO^ and lAbbasi et al.l toilh . 
Also shown are sen s itiviti es (so lid and dashed lines) for the current analy sis and for the search 
from lAbbasi et al.l (120111 ) (see lAbbasi et al.l (120091 ): iThrane et al.l (120091 ) for the sensitivities 
corresponding to those searches). 



The present limits are more stringent than those obtained for the northern hemisphere 
by previous multi year experiments (also indicated in the figure) and are competitive with 
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Table 1: Results of the candidate list search. The source coordinates and the Q and p- 
values are shown as well as the limits on the flux intensity ^j^^ latter has units 

10~^GeV~^cm~^s~^. The rows are sorted in order of increasing p-value. 
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those set by the IceCube collaboration ( lAbbasi et al.l 1201 ll ) for 6 < —30°. It should be 
noted that even though they all set limits on the intensity of an assumed spectrum, the 
experiments are sensitive in different energy ranges. For this spectrum, ANTARES detects 
most events at energies in a broad range around 10 TeV, which is the relevant energy range 
for several galactic sources (jCrocker et al.ll2005l ). Northern hemisphere limits shown from 
the IceCube experi ment probe the neutrino flux predominantly in the region above 1 PeV 
jAbbasi et al.ll201l[ l. 

The event selection and the search method have been cross-checked with an independent 
analysis using the same s election criteria and a search method based on the ex pectation- 
maximization algorithm (jPempster et al.l 119771 : lAguilar &: Hernandez- Reyl l2008l ). In this 
method, the angular spread of the signal events is a free parameter in the likelihood and the 
maximization is performed analytically. The results of both the full sky and the candidate 
list search are consistent with the results discussed earlier. 



6. Autocorrelation Analysis 

In a separate analysis, the cumulative distribution of the number of event pairs as a 
function of their angular separation has been studied using the neutrino candidate events. 
This autocorrelation analysis has been employed to search for features in the data such as 
an excess of tracks at an arbitrary level of angular separation. The study is independent of 
detector simulations and neutrino source models. 

A reference autocorrelation distribution is determined by scrambling the data itself 
approximately one million times and averaging the resulting distributions. The comparison 
between the data and the reference distribution is performed following Li and Ma (1983). 

The maximum excess above the reference distribution has a significance of 1.1 a at 
angular scales smaller than 7 degrees. Such a deviation is expected from a collection of 
random background events with a (trial factor corrected) probability of 55%. 



7. Conclusions 

A search for cosmic sources of high energy neutrinos using the data taken with the 
ANTARES neutrino telescope during the first two years of operation has been presented. 
For the first half of the data considered, the detector was about half its final size. An estimate 
of the angular resolution constrained by data yields 0.5 ± 0.1 degrees (0.4 degrees for the 
12-line detector), confirming the expectation that excellent angular resolution is achievable 
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using sea water as the detection medium. Neither the full sky search nor the candidate list 
search show a significant excess of events. Limits have been obtained on the high energy 
neutrino fiux for a number of selected source candidates. For many candidate sources, the 
limits presented here are the most stringent to date. 
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